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Abstract 

Recently, some efforts focus on differentiating dark energy and modified gravity with the growth func- 
tion S{z). In the literature, it is useful to parameterize the growth rate / = dlnS/dlna = flZi with 
the growth index 7. In this note, we consider the general DGP model with any Qk- We confront the 
growth index of DGP model with currently available growth rate data and find that the DGP model is 
still consistent with it. This implies that more and better growth rate data are required to distinguish 
between dark energy and modified gravity. 
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1 Introduction 

The current accelerated expansion of our universe [Tl[2l[3l[4l[5l[6l[7l[8l[9l |49] has been one of the most 
active fields in modern cosmology. There are very strong model-independent evidences [TO] (see also e.g. [TT| ) 
for the accelerated expansion. Many cosmological models have been proposed to interpret this mysterious 
phenomenon, see e.g. [T] for comprehensive reviews. 

In the flood of various cosmological models, one of the most important tasks is to discriminate between 
them. Recently, some efforts have been made. For instance, it is important to determine that the dark 
energy is cosmological constant or dynamical dark energy [I]. Caldwell and Linder proposed a so-called 
w — w' analysis in 12J to discriminate dark energy models, and then was extended in [131 114j . A recent 
review on w — w' analysis can be found in |15j . Another tool to discriminate models is the statefinder 
diagnostic proposed by Starobinsky et al. in [16j . For a comprehensive list of relevant works on w — w' 
analysis and statefinder diagnostic, one can see e.g. [TJ and references therein. 

Recently, some efforts to discriminate models focus on differentiating dark energy and modified gravity 
with the growth function 5[z) = Spm/ Pm of the linear matter density contrast as a function of redshift z. By 
now, most of cosmological observations merely probe the expansion history of our universe [U [21 El HI El El 
El El El |49j . As well-known, it is very easy to build models which share a same cosmic expansion history by 
means of reconstruction between models. Therefore, to distinguish various models, some independent and 
complementary probes are required. Recently, it is argued that the measurement of growth function 6{z) 
might be competent, see e.g. [IHlIiaEniEIlEaEaiaESlEaEaEHlEaEniEIlEHI-Iftwo models, especially 
dark energy and modified gravity models, share a same cosmic expansion history, they might have different 
growth history. Thus, they could be distinguished from each other. 

One of the leading modified gravity models is the so-called Dvali-Gabadadze-Porrati (DGP) braneworld 
model [321 ES] J which altering the Einstein-Hilbert action by a term arising from large extra dimensions. For 
a list of references on DGP model, see e.g. [211 [31] and references therein. The first approach to study the 
growth function 5{z) of DGP model is numerical solution, see e.g. [23 [27l [35] . The second approach is to 
parameterize the growth rate / = dlnJ/dlna, where a = (1 -f z)~^ is the scale factor of our universe. 

For many years, it has been known that a good approximation to the growth rate /, within Einstein 
gravity, is given by 

dm a 

where 7 is the growth index, whereas is the fractional energy density of matter. In the very beginning, 
Eq. ([T]) was introduced in [361137] . There, / was defined purely in terms of the present value, using the present 
matter density, not valid for arbitrary redshift. Also, it was not until [21j that it was applied to anything 
beyond matter, curvature, and a cosmological constant. Finally, not until [21] was it applied to gravity 
other than general relativity, and then in [TH] generalized to modified gravity, varying equation of state, and 
an integral relation for growth. This parameterized approach has been tested in some works recently, see 
e.g. [TSl[Il[2ni[21[2i[2Sll3nil3Ill3E]- The theoretical value of 7 for ACDM model is 6/11 ~ 0.545 [H [H] 
whereas 7 ~ 0.55 for other parameterized dark energy models [H]. The theoretical growth index of flat DGP 
model (whose fife = exactly) is 7 = 11/16 = 0.6875 [H] (In fact, 7 = 11/16 is a high redshift asymptotic 
value. 7 = 0.68 is favored for the fit to the whole growth history to the present, while another approximation 
is given in Eq. (27) of [15] for any redshift, with 7 — 7/11 in the asymptotic future). Therefore, it is possible 
to distinguish the dark energy model (including ACDM model) from the flat DGP model. 

In this note, we consider the general DGP model with any Vlk- Wc confront the growth index of DGP 
model with currently available growth rate data and find that the DGP model is still consistent with it. On 
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the other hand, it is shown that the ACDM model is also consistent with current growth rate data [21] ■ This 
implies that more and better growth rate data are required to distinguish between dark energy and modified 
gravity models. 



2 Growth index of DGP model 

In the DGP model, the Friedmann equation is modified as [33] (see also e.g. [28 ] \29 \ l39]) 



Tr V 3 



Pirn 



(2) 



where H = a/a is the Hubble parameter; a dot denotes the derivative with respect to cosmic time t] the 
constant j'c is the crossover scale; spatial curvature fc = 0, fc > and fc < for flat, closed and open universe, 
respectively. Here, we only consider the self-accelerating branch. It is easy to rewrite Eq. ^ as an equivalent 
form [33] (see also e.g. [iU] ) 

(3) 
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where pr^ = 3/(327rG'r^). One can alternatively recast Eq. |2]) as the "standard" form 
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in which we have used Eq. ([3]) in the last equality. Here, pde can be regarded as an effective "dark energy" 
component, which compiles the contributions to the Friedmann equation from the extra dimensions |281 129j . 
From Eq. we have 

1 - r^fc = fl,n + ^de, (6) 

where flk = —k/{a'^H^), and fli = {8'!TGpi)/{3H^) for i — m and de. From Eqs. (O, ([6]) and the energy 
conservation equation p^e + 3H{pde + Pde) = 0, we find that the effective equation-of-state parameter (EoS) 
for the effective "dark energy" component Wde = Pde/ Pde is given by 



Wde 



(7) 



For the flat DGP model whose ilk — exactly, we get Wde — ^(1 + ^m) ^, which agrees with the one 
of [281 [29]. From Eq. ([3]), it is easy to see that 



H 



(8) 



where we have used pm = Pmoo. the subscript "0" indicates the present value of the corresponding quantity; 
Qr, = (SnGprJ/iSH^) = l/iAH^rl). From Eq. we have 



1 = 



kO- 



(9) 



Therefore, only two of Qmo, ^-nd ilko are independent model parameters. For the flat DGP model whose 
i^fco = exactly, Eq. @ becomes f^„iO — ^ ~ 2-v/iX~. 
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In Einstein gravity, the growth function 6{z) at scales much smaUer than the Hubble radius obeys the 
following differential equation [HI HSl EOl EH [HI SO] 

6 + 2HS = AnGpmS. (10) 

In modified gravity, Eq. ^ has been modified to [271 US ES HO] 

5 + 2H5 = ATTG,sPmS, (11) 

where GeS is the effective local gravitational "constant" measured by Cavendish- type experiment, which is 
time-dependent. In general, GeS can be written as 

Ge. = G(l+i^), (12) 

where f3 is determined once we specify the modified gravity theory. In the DGP gravity, /3 is given by [27j 

ESKSaSO] 

P=l-2r,H{l + ^]. (13) 



By using Eq. ([3]), /3 can be written as [40] 



3/3 sni-sii-rikf + a' ' 



where 



a = 2v/i - r^fc (3 - 4f]fe + 2n^nk + nl) . (15) 

For the flat DGP model whose fik — exactly, we get /3 = —(1 + il^„)/(l — ilf„), which agrees with the one 
of |28lE9]. Following [2Tl|3l], we rewrite Eq. HH) as 

(In^)" -t- {ln6f + (2 + 1^) (ln<5)' = ^ (l + ^) ^m, (16) 

where a prime denotes the derivative with respect to In a. By using Eqs. ^ and the energy conservation 
equation pde + 3H{1 + Wde)pde — 0, we find that 

f = ;^ = -^ + ^-^^.e(l-0,-an), (17) 



where Wde is given in Eq. ([7]). Substituting Eq. ([T7)) into Eq. p^ . we obtain 
(ln<5)" + (ln(5)" + (ln5)' 



1 3 

- (1 + ilfe) - -Wde [l-^k~ ^m) 



+ (18) 



In fact, the growth rate / = dln^/dlna = (In (5)'. By using the definition of fJ,,,,, the energy conservation 
equation + iHpm = and Eq. (fTT]) , we have 

^'^ = [iWde (1 - f^fc - f^m) " ^^fc] ■ (19) 

Therefore, we find that 

d f 

{\n5)" ^f'^flra [3Wde {I ~ ^^k ^ ^m) - ^k] . (20) 
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So, Eq. (fTH]) becomes 

[3wde (1 - rife - - ^k] -r^ + f + f 



1 3 

- (1 + ilfe) - -lUrfe (1 - f^fc - fir. 



^ ^ + -5- ) (21) 



2 V 3/3 

where Wde and are given in Eqs. ([7]) and p^ . respectively. Substituting Eq. (IT|) into Eq. (PT|) and expanding 
around ri,„ 1 (good approximation especially at z>l), after some tedious algebra, we finally arrive at 

107 {n„^ - 1) - rife (1 - 27) + 3 (1 - 27) (1 - a„) = 6 - 1) + - 1) (a„ + 1) + 2nk] , (22) 

where we have ignored the higher order terms of small quantities 1 — flm and 17^.. Noting that Eq. namely 
1 — flm = flk + fide, wc cousidcr three cases for Eq. (|22p . 



Case (I) ilfc = exactly or fife ^ fi^e = 1 — ftm <C 1. 

In this case, throwing out the terms of flk in Eq. ([22|) . then eliminating (fi™ — 1) in both sides, and 
using ftm 1 finally, we have 7 = 11/16. Obviously, it agrees with the known one of the flat DGP 
model [TBI. 

Case (II) fide < fife = 1 - fim < 1. 

In this case, eliminating {flm ^ 1) = ^flk in both sides of Eq. (|22p and then using flm — 1, we have 
7 = 4/7. In fact, it coincides with the curvature solution found firstly in |51] and could be considered 
as a special case of the results in [2TJ [18] with w — —1/3. 

Case (III) fife ^ fide ^ I - flm < 1- 

In this case, fife, fi^e and 1 — fi^, are at the same order. Noting that 1 — fi„i = fife + fi^e, for 
convenience, we parameterize fife = m(l — fim), where < to < 1. It is worth noting that generally 
TO is time-dependent and one considers only an instantaneous value of m. Substituting into Eq. (|22p. 
then eliminating (fi^ — 1) in both sides, and using flm — > 1 finally, we have 

11 — 3to 

Obviously, when to — and 1, 7 — *■ 11/16 and 4/7, respectively. 

For the fiat DGP model whose fife = always, the only theoretical growth index is given by 11/16. For 
the DGP models whose fife 7^ 0, the situations are different. Noting that fim c>c (1 + z)'^, fife c>c (1 + z)^ 
and fide oc (1 + 2;)3(i+™de) -vvith Wde < —1/3 to accelerate the expansion of our universe, fife increases faster 
than fide when z increases. So, for high z, fi^e ^ fife = 1 — fim ^ 1 eventually. Thus, we have 7 — 4/7 
eventually, regardless of the value of fifeo- However, if fifeo deviates from very small, it is still possible that 
fife <C fide = 1 — flm ^ 1 at high z, where 7 = 11/16; eventually 7 = 4/7 at higher z. On the other hand, if 
fifeo deviates from not so small, the only theoretical growth index is given by 4/7. 

Finally, it is worth noting that the above results are obtained as high redshift asymptotic values, rather 
than values for the whole growth history or the asymptotic future. On the other hand, the general solution 
of 7 is obtained in [18] for the fiat DGP model. We refer to [18] for details on this. 



3 Confronting with current growth rate data 

The most useful currently available growth rate data involve the redshift distortion parameter observed 
through the anisotropic pattern of galactic redshifts on cluster scales [12] (see also |3T]). The parameter Pl 
is related to the growth rate / as [31] 

Pl = {, (24) 
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z 


(iL 


h 


fobs 


Reference 


0.15 


0.49 ±0.09 


1.04±0.11 


0.51 ±0.11 


m 


0.35 


0.31 ±0.04 


2.25 ±0.08 


0.70 ±0.18 


m 


0.55 


0.45 ±0.05 


1.66 ±0.35 


0.75 ±0.18 


m 


0.77 


0.70 ±0.26 


1.3±0.1 


0.91 ±0.36 


m 


1.4 


o.eot^jt 


1.5 ±0.20 


0.90 ±0.24 


m 


3.0 






1.46 ±0.29 


m 



Table 1: The currently available data for and b at various redshifts, along with the inferred growth rates. 
Notice that Ref. [47] only reports the growth rate and not the (3l and b parameters, since the growth rate was 
obtained directly from the change of power spectrum Ly — a forest data in SDSS at various redshift slices. 
This is an extended version of the dataset used in Ref. [31], and contains a new data point at z = 0.77 |45j . 



where 6 is the bias factor. We present the currently available data for (3l and b at various redshifts in Table[l] 
along with the inferred growth rates. This is an extended version of the dataset used in Ref. [31 , and contains 
a new data point at z — 0.77 [45]. Hereafter, we call the five data points except the one at z = 0.77 as dataset 
"Fobs" and call all the six data points as dataset "Fobsext" . 

About the currently available growth rate data, it is worth noting that there is considerable variation in 
analysis of different references and hence there is no accomplished consensus in fact. On the other hand, 
many of these results are not of the growth rate in isolation but fold in other large scale structure information. 
Besides, in Table [1] some of the quoted values include the other datasets in the analysis, e.g. [46j includes 
[i5] to derive their value of So, one should use these data with caution. 

Using the growth rate data in Table (T] we can perform a analysis to find the growth index 7 and check 
its consistency with the theoretical values. As well-known, the corresponding reads 

2 



X^(P,7)=E 



fobsiZt) - fthiz^;P,J) 



(25) 



where fobs and its corresponding Icr uncertainty Uf^^^ are given in Table [U p denotes the model parameters; 
fth{zi', p,7) can be obtained from Eq. ([1]), in which Qm can be rewritten as a more convenient form 

n..M - ^^i^. (26) 

and E{z) can be found in Eq. ([5]). 

It is worth noting that the references in Table [T] have assumed flat ACDM (with ilmO.A = 0.30 for the 
five data points except the one at z = 0.77, and with ilmO,A = 0.25 for the data point at z = 0.77) when 
converting redshifts to distances for the power spectra and therefore their use to test models different from 
ACDM might be unreliable, as stressed in [31]. However, we can get around this problem in a new way. As 
mentioned above, the key is the redshift-distance relation. If the DGP model and ACDM model share the 
same redshift-distance relation, these growth rate data can also be used in the DGP model. To this end, we 
should properly select the model parameters of DGP model in order to reproduce the same redshift-distance 
relation of ACDM model. There is a simple and efficient method. As well-known in any textbook, the 
comoving distance is given by 

r(z) ^ — ^ F ( r J^) , (27) 



HoVWi V Jo m) 
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where the function F(x) = x, sinx and sinhx for 51^ = 0, fit < and f2fc > 0, respectively. For convenience, 
we use the dimensionless comoving distance Hor{z) instead. It is easy to get the Hor{z) hnc for the flat 
ACDM model with flmO,A = 0.30. Then, we discretize it into many points in a redshift range (for instance, 
< z < 4, which covers the range of current growth data), and manually assign a relative "error" (say, 
0.5%) to these discrete points. So, we have many fake "data points" in hand. Then, we fit the DGP model 
to these fake "data points" and find out the "best fit" parameters. Obviously, the DGP model with these 
"best fit" parameters will share the same redshift-distance relation with the flat ACDM model. Hence, the 
growth rate data in Table [1] can also be used for the corresponding DGP model. Finally, there is a minor 
remark on the procedure mentioned above. In fact, it holds only for the low redshift quantities considered 
here. Comparison of the distances to CMB last scattering (z ~ 1100) will show large differences between 
these matching models. Fortunately, the redshift of current growth data is less than 4 and the procedure 
mentioned above works well. 




0.5 1 1.5 2 2.5 3 3.5 4 0.5 1 1.5 2 2.5 3 3.5 4 

z z 

Figure 1: The Hor{z) of flat DGP model with the corresponding "best fit" parameters (red dashed line) and 
flat ACDM model with r2mO,A — 0.30 (black solid line) are shown in left panel. The right panel shows the 
relative departure between these two Hgr^z) lines. 



We consider the flat DGP model whose flka = exactly at flrst. Following above method, we flnd that 
the "best fit" parameter of flat DGP model is fi^c — 0.1519 (the corresponding ri„io — 0.2204). In Fig. [l] 
we show the Hor{z) of the corresponding flat DGP model along with the one of flat ACDM model with 
^mO.A — 0.30. The largest relative departure between these two HQr{z) lines in range < z < 4 is less than 
1.5%. Similarly, we consider the general DGP model with any flko and flnd that the "best flt" parameters 
of DGP model are n^o = 0.2217 and fi^^ = 0.1734 (the corresponding ilko = -0.0921). In Fig. [21 we show 
the Hor{z) of the corresponding DGP model along with the one of flat ACDM model with flmO,A = 0.30. 
Obviously, these two Hor{z) lines are degenerate in fact. The largest relative departure between the Hor{z) 
of DGP model with the corresponding "best fit" parameters and fiat ACDM model with fimO.A — 0.30 in 
range < z < 4 is less than 0.19% surprisingly. The DGP model with corresponding "best fit" parameters 
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excellently reproduces the same redshift-distance relation of the flat ACDM model with r^mO.A — 0.30. So, 
the growth rate data in Table [1] can also be used for the corresponding DGP model. It is worth noting that 
all the above "best fit" parameters are consistent with the constraints of [2211151 ISH] (see also e.g. [5H [5U] ) 
and therefore they are also observationally acceptable. 




Here, we fit the growth index 7 of these "reproduced" DGP models to the growth rate data in Table [1] 
For the case of flat DGP model whose O^o = exactly, setting Slmo and fJ^e as the corresponding "best 
fit" values and minimizing ^ with respect to 7, we find that for dataset Fobs, 7 — 0.438^0 ^i^^ with 1<t 
uncertainty, or 7 = G.438_q 209 with 2a uncertainty; for dataset Fobsext, 7 — 0.429_Q]^Qg with Icr uncertainty, 
or 7 = 0.4291q'2Q5 with 2a uncertainty. Obviously, the only theoretical growth index 7 = 11/16 = 0.6875 is 
consistent with these results at 2a uncertainty. 

For the case of DGP model with fifco 7^ 0, setting ilmOj and fifco as the corresponding "best fit" values 
and minimizing with respect to 7, we find that for dataset Fobs, 7 — QA%h'^Q^{'^ with Icr uncertainty, 
or 7 = 0.465^0 221 with 2a uncertainty; for dataset Fobsext, 7 = 0.4571° iji^g with \a uncertainty, or 7 = 
0.457lo'2i7 with 2a uncertainty. Obviously, the theoretical growth index 7 — 11/16 = 0.6875 is consistent 
with these results at 2a uncertainty, whereas the other theoretical growth index 7 = 4/7 ~ 0.5714 is consistent 
with these results at la uncertainty. 



4 Concluding remarks 

In summary, we consider the growth index of general DGP model with any Qk in this note. We confront the 
growth index of DGP model with current growth rate data and find that the DGP model is still consistent 
with it. On the other hand, it is shown that the ACDM model is also consistent with current growth rate 
data [21] . This implies that more and better growth rate data are required to distinguish between dark energy 
and modified gravity models. 



9 



Some remarks are in order. Firstly, it is worth noting that the references in Tabled] have assumed flat 
ACDM (with r^mO.A = 0.30 for the five data points except the one at z = 0.77, and with r^mO.A = 0.25 for the 
data point at z = 0.77) when converting redshifts to distances for the power spectra and therefore their use 
to test models different from ACDM might be unreliable, as stressed in |3T]. Although we can get around this 
problem by means of reproducing the same redshift-distance relation between the DGP model and ACDM 
model, the really complete solution to this problem lies on the reanalyzing the power spectra with the proper 
redshift-distance relation in the corresponding modified gravity. Secondly, in fact, we fixed Qmo while fitting 
for 7 in section [31 This is required to ensure the validity of matching procedure of finding an equivalent 
DGP model. We admit that this weakens our conclusion. It is desirable to find a new method to improve 
the validity of matching procedure in the future works. Thirdly, one may use other parameterizations to the 
growth rate /, such a,s f — ri7n(l + v) 00], or / = il^ with 7 = 7o + JqZ [38\. However, the additional 
parameters and 7q are found to be negligible in fact [SD] [SH] ■ On the other hand, since currently available 
growth rate data points are so few and have large errors, it is difficult to tightly constrain these additional 
parameters. Fourthly, although as shown in [23j that non-trivial dark energy clustering or interaction between 
dark energy and dark matter might bring some troubles, we consider that combining the probes of expansion 
history and growth history is still promising to distinguish the dark energy and modified gravity. Of course, 
new idea to distinguish dark energy and modified gravity is still desirable. 
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